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Abstract The CCD photometric observations of the eclipsing binary PS Persei (PS Per) 
were obtained on two consecutive days in 2009. The 2003 version Wilson-Devinney code 
was used to analyze the first complete light curves in V and R bands. It is found that PS 
Per is a short-period Algol-type binary with the less massive component accurately filling 
its inner critical Roche lobe. The mass ratio of q = 0.518 and the orbital inclination of 
i = 89.°86 are obtained. On the other hand, based on all available times of primary light 
minimum including two new ones, the orbital period has been improved. 
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1 INTRODUCTION 

PS Per (a 20 oo.o = 02 h 39 m 33. ;5 3 and £2000.0 = +45°38'05."5) was designated by Kukarkin et al. 
(1968). But photographic and visual times of light minimum have been obtained since 1926. Later, 
Photoelectric and CCD times of light minimum were published by Safaf & Zejda (2000a), Safaf & 
Zejda (2000b), Zejda (2002), Agerer & Hubscher (2003), Zejda (2004), Diethelm (2005), Hiibscher et 
al. (2006), Zejda et al. (2006), Brat et al. (2009), and Diethelm (2010). But, no complete light curve of 
the binary system have been made so far for photometric analysis. 

In this paper, the first complete light curves in V and R bands were presented. And the absolute 
physical parameters as well as orbital period were determined. 

2 OBSERVATIONS 

New CCD photometric observations of PS Per in V and R bands were carried out on 2009 November 
13 and 14 using the 85-cm telescope at the Xinglong Station of National Astronomical Observatory of 
China (NAOC), equipped with a primary-focus multicolor CCD photometer. The telescope provides a 
field of view of about 16. 5 x 16. 5 at a scale of 0. 96 per pixel and a limit magnitude of about 17 mag 
in V band. 

The typical exposure times in V and R bands were 90s and 60s respectively. The coordinates of the 
variable, comparison, and check stars are listed in Table 1 . The data reduction was performed by using 
the aperture photometry package IRAfQ (bias subtraction, flat-field division). Extinction corrections 
were ignored as the comparison star is very close to the variable. In total, 445 CCD images in the V 
band and 446 images in the R band were obtained. Several new times of light minimum (see Table 2) 
are derived from the new observation by using a parabolic fitting method. 

The first complete light curves in V and R bands are obtained, and displayed in the top panel of 
Figure 1. The new orbital period revised in the next section was used to calculate the phase. 



IRAF is developed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under contract to the National Science Foundation. 
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Table 1 


Coordinates of PS Per and its Comparison and Check Stars. 






Stars 


Q2000 


52000 






PS Per 


02 ,l 39 m 33.3 s 


45°38'05.5" 






Comparison 


02 ,l 39 m 24.1 s 


45°42'22.1" 






Check 


02' t 39 m 22.2 s 


45°43'34.0" 




Table 2 


New CCD Times of Light Minimum for PS Per. 




No. 


J.D. (Hel.) (days) 


Error (days) 


Min. 


Filter 


1 


2455149.2742 


±0.0004 




V 




2455149.2752 


±0.0005 




R 


2 


2455149.9768 


±0.0003 




V 




2455149.9768 


±0.0004 




R 


3 


2455150.3267 


±0.0004 


II 


V 




2455150.3265 


±0.0004 


II 


R 



0.0 






0.0 


0.2 


0.4 0.6 

Phase 


0.8 


1.0 
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Fig. 1 Top panel: the light curves of PS Per in the V 
and R bands obtained on 2009 November 13 and 14. 
The points in R band has been shifted down by 1.0 
mag. Bottom panel: the differential light curves of the 
comparison star relative to the check star. The points in 
R band has been shifted up by 0.2 mag. 

3 ORBITAL PERIOD STUDY 

All available times of primary light minimum seen in literature were collected and listed in Table 3, 
which also includes the data in the database of eclipsing binaries of Kreiner (2004). For my two-band 
light minima, a mean time of light minimum is given. The O — C values of all minimum times were 
computed with the ephemeris given by Kreiner et al. (2001): 



Min I = 2424527.2165 + d . 70217968 x E, 



(1) 
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and listed in the fifth column of Table 3. The corresponding O — C diagram, Figure 2, shows that the 
data are distributed around a straight line. So, a linear ephemeris was used to fit the O — C values. The 
photographic and visual data show large deviation from the straight line for their low quality. In the 
fitting process, a weight of 10 is used for photoelectric and CCD data, and 1 for photographic and visual 
data. The CCD times, 2451841.3121 and 2452213.4662, have the weight of 1 for their large errors. A 
least-squares fit to the data gave the following ephemeris: 

Min I = 2424527.2163 + d . 70217977 x E, (2) 

The new ephemeris is plotted in Figure 2 with a solid line. The residuals with respect to Equation (2) 
are listed in the sixth column of Table 3. 
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Fig. 2 O — C diagrams. The solid line is calculated 
with the new ephemeris in Equation (2). 



4 PHOTOMETRIC SOLUTIONS WITH THE W-D METHOD 

The light curves were analyzed using the 2003 version of the Wilson-Devinney code (Wilson & 
Devinney 1971; Wilson 1979, 1990). Since the spectral type of PS Per is F5, an effective tempera- 
ture of T\ = 6750K is assumed for the primary component. Assuming the photospheric surface of 
the binary star is convective, gravity-darkening coefficients (171 = 92 = 0.320) and bolometric albedo 
(A\ = A2 = 0.5) were used. According to the tables of van Hamme (1993), the limb-darkening coeffi- 
cients 0.506 for V band (xi V = 0.506) and 0.414 for R band (xi R = 0.414) were adopted. 

Since no mass ratio has been published in literature, a g-search method was used to determine the 
mass ratio. Solutions were carried out for a series of values of the mass ratio q — M2/M1 (q = 0.3, 0.4, 
0.5, 0.6, 0.7, 0.8, 0.9, 1.0). Considering the light curves of EB type, mode 2 (detached mode) is assumed. 
The behavior of the sum of the residuals squared, S, as a function of mass ratio q is plotted in Figure 
3, showing that £ reaches the minimum value near q = 0.5. Therefore, the mass ratio was taken as 
an adjustable parameter and given the initial value of q = 0.5. After some differential corrections, The 
solution converged to mode 5 (semi-detached) and gave the final mass ratio of q = 0.518. The derived 
physical parameters are listed in Table 4. The theoretical light curves computed with the parameters are 
plotted in Figure 4 as a solid line. 
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Fig. 3 Relation between £ (the sum of the residuals 
squared) and q for PS Per. 
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Fig. 4 Same as the top panel of Figure 1 . But the solid 
curves represent the theoretical light curves computed 
with the parameters in Table 4. 

5 DISCUSSION AND CONCLUSIONS 

In this paper, my photometric solution reveals that PS Per is a semi-detached system. The Roche- 
geometry configuration that the less massive and cool secondary component fills its inner Roche lobe 
permits a dynamical mass transfer from the secondary to the more massive primary star, suggesting a 
continuous period increase just as in AI Cru (Zhao et al. 2010) and DD Mon (Qian et al. 2009). The 
orbital period of PS Per, however, does not show continuous increase in this paper. This may be due to 
the low quality of photographic and visual times, and the short span of photoelectric and CCD times. In 
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order to confirm the mass transfer from the secondary to the primary star, long-term orbital timing data 
are required. 

According to the derived physical parameters listed in Table 4 and the Harmanec's (1988) relation 
for masses and radii as functions of spectral type, the following orbital parameters can be derived: 
M x = 1.31 Mq, Rx = 1.39 Rq, Li = 3.62 L Q , M 2 = 0.68 M Q , R 2 = 1.35 R Q , L 2 = 0.89 L Q , and 
a = 4.19 Rq. In order to further verify the parameters, Spectroscopic observations of the radial velocity 
curves of both components are needed. 

As showed in Figure 1, the second maxima of the light curves are a little higher than the primary 
maxima. The weak O'Connell effect maybe arises from a hot spot on the primary component as a result 
of the impact of the gaseous stream from the cooler, less massive secondary component. Such hot spot is 
often seen in other semi-detached binary systems, such as CL Aur (Lee et al. 2010) and KQ Gem (Zhang 
2010). Considering the late spectral type and fast rotation of the secondary star, The asymmetry of the 
light curves can be also attributed to a cool spot on the secondary star caused by magnetic activity. It is a 
reasonable trial that the magnetic activity makes the light curves show more variability than the impact 
of the gaseous stream does. So, in order to tell the two mechanisms of magnetic activity and impact of 
the gaseous stream, the investigation on long-term behaviour of the light curves is also needed. 

Acknowledgements I thank an anonymous referee for some useful suggestions. This work is supported 
by Natural Science Foundation of Shanxi Normal University (No. ZR09002). 

References 

Agerer, E, & Hubscher, J. 2002, IBVS, No. 5296 
Agerer, E, & Hubscher, J. 2003, IBVS, No. 5484 
Brat, L., Trnka, J„ & Lehky, M. 2009, OEJV, No. 107 
Diethelm, R. 2003, IBVS, No. 5438 
Diethelm, R. 2004, IBVS, No. 5543 
Diethelm, R. 2005, IBVS, No. 5653 

Hubscher, J., Paschke, A., & Walter, F 2006, IBVS, No. 5731 
Diethelm, R. 2010, IBVS, No. 5920 

Harmanec, P. 1988, Bull. Astron. Inst. Czechoslovakia, 39, 329 

Kukarkin, B. V, Efremov, Yu. N., Frolov, M. S., et al. 1968, IBVS, No. 31 1 

Kreiner, J. M., Kim, C.-H, & Nha, I.-S. 2001, An Atlas of O — C Diagrams of Eclipsing Binary Stars (Krakow: 

Wydawnictwo Naukowe Akad. Pedagogicznej) 
Kreiner, J. M. 2004, Acta Astron., 54, 207 

Lee, J. W., Kim, C.-H., Kim, D. H, Kim, S.-L., Lee, C.-U., & Koch, R. H. 2010, Al, 139, 2669 

Qian, S.-B., Zhu, L.-Y, Boonrucksar, S., Xiang, F.-Y, & He, J. -J. 2009, PASI, 61, 333 

Safaf, J., & Zejda M. 2000a, IBVS, No. 4887 

Safaf, J., & Zejda M. 2000b, IBVS, No. 4888 

van Hamme, W. 1993, A J, 106, 2096 

Wilson, R. E., & Devinney E. J. 1971, Apl, 166, 605 

Wilson, R. E. 1979, Apl, 234, 1054 

Wilson, R. E. 1990, Apl, 356, 613 

Zejda, M. 2002, IBVS, No. 5287 

Zejda, M. 2004, IBVS, No. 5583 

Zejda, M., Mikulasek, Z., & Wolf, M. 2006, IBVS, No. 5741 
Zhang, L.-Y. 2010, PASP, 122, 309 

Zhao, E.-G., Qian, S.-B., Fernandez Lajus, E., von Essen, C, & Zhu, L.-Y. 2010, RAA(Research Astron. 
Astrophys.), 10, 438 



6 



Jinzhao Yuan 



Table 3 Times of primary light Minimum for PS Per. 



JD (Hel.) 


Method 


Error 


E 


O — C 


Residuals 


Ref. 


2424527.250 


P 







0.03350 


0.03374 


MVS 2 


2428834.407 


P 




6134 


0.02034 


0.02000 


MVS 2 


2430972.531 


P 




9179 


0.00722 


0.00660 


MVS 2 


2430991.464 


P 




9206 


-0.01863 


-0.01925 


MVS 2 


2435718.545 


P 




15938 


-0.01124 


-0.01249 


MVS 2 


2435725.534 


P 




15948 


-0.04404 


-0.04529 


MVS 2 


2436114.558 


P 




16502 


-0.02758 


-0.02888 


MVS 2 


2436596.315 


P 




17188 


0.03416 


0.03278 


MVS 2 


2436603.330 


P 




17198 


0.02736 


0.02598 


MVS 2 


2436850.460 


P 




17550 


-0.00988 


-0.01128 


MVS 2 


2436852.580 


P 




17553 


0.00358 


0.00217 


MVS 2 


2436876.438 


P 




17587 


-0.01253 


-0.01393 


MVS 2 


2436895.406 


P 




17614 


-0.00338 


-0.00479 


MVS 2 


2437588.460 


P 




18601 


-0.00073 


-0.00223 


MVS 2 


2437939.528 


P 




19101 


-0.02257 


-0.02412 


MVS 2 


2437944.477 


P 




19108 


0.01117 


0.00961 


MVS 2 


2437946.551 


P 




19111 


-0.02136 


-0.02291 


MVS 2 


2437970.442 


P 




19145 


-0.00447 


-0.00602 


MVS 2 


2438321.547 


P 




19645 


0.01069 


0.00908 


MVS 2 


2438385.453 


P 




19736 


0.01834 


0.01672 


MVS 2 


2447028.565 


V 




32045 


0.00065 


-0.00211 


BBSAG 85 


2447118.446 


V 




32173 


0.00266 


-0.00011 


BBSAG 86 


2447170.406 


V 




32247 


0.00136 


-0.00142 


BBSAG 87 


2447384.566 


V 




32552 


-0.00344 


-0.00625 


BBSAG 89 




V 
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dL /U4 


-U.UUO / J 


-U.UUbO / 


DDCAr on 

rSrJoALr y\) 


7J.A7^AA All 
Z44 / jOO-4jZ 


V 




"37Q 1 1 
jZo 1 1 


n 001 qs 
-u.uuiyo 


O OO/IS 1 
-U.UU45 1 


Aw y 1 


Z44 / o!?4. j4o 


V 




DDL 1 


o oo^sq 


O OOA77 
-U.UUO/ / 


ddc An QA 
ljo/Wj yH- 


DAASZ 1 ^Q7 
Z44ol JO.J7 / 


V 




DDOZ.D 


o ooass 

-U.UUO00 


O OOQ7Q 
-U.UUV ly 


ddc An OA 


7/1/1 Q7Q3 3A7 
Z44oZOj. jOZ 


V 


.UU4 


ddodZ 


U.UUZO / 


O OOOTA 

-U.uuuoo 


ddc An Q7 


7/1/lG^OQ A^Sl 
Z44ojUy.40o 


V 


.UUJ 


1A\ ZA 
041 J4 


n oot>q 




ddc An QQ 


ZH-'+ooO / . J / D 


V 


.UUJ 


lAfsfxA 
D^tVOH 


o nnnn7 

u.uuuu / 


O 007QA 

-u.uuzyn- 


RRC An 1 0? 

DDonU 1UZ 


7/1AQ707 ^ 1 1 
Z44VZUZ. J 1 1 


V 




^ 1 A\ 
DD 141 


noi 

-U.UUIOj 


O OOzlAfi 
-U.UU4DO 


ddc An 1 C\A 
DDonU 1U4 


lAAQ^Afi ^7^ 
Z44VJ40.J / J 


V 


.uuo 


1 

jjOj 1 


-U.UUJOo 


O OOS7S 
-U.UUO / 


ddc An 1 07 

ddjAu 1U / 


Z44VDJj.j1Z 


V 




DD 1 oD 


00001 
u.uuuu 1 


00^ 1 
-U.UUj 1U 


ddc An 1 os 

DDonU lUo 


lAAQQA^ AC\\ 
Z44VV4 J .4U 1 


V 


.UUJ 


Q^l QQ 


O OI 11 A 
-U.U1 / /4 


O 070SQ 
-U.UZU07 


ddc An 1 1 
rSrSoAO 1 1U 


2450713.6127 


ec 


.0014 


37293 


0.00939 


0.00612 


Safar & Zejda (2000a) 


2450721.3309 


cc 


.0008 


37304 


0.00362 


0.00035 


BBSAG 116 


2450839.2981 


cc 


.0021 


37472 


0.00463 


0.00135 


Safar & Zejda (2000b) 


2450841.4035 


cc 


.0021 


37475 


0.00349 


0.00021 


Safar & Zejda (2000b) 


2451077.3366 


cc 


.0006 


37811 


0.00422 


0.00091 


BBSAG 119 


2451088.572 


V 


.003 


37827 


0.00474 


0.00142 


BBSAG 119 


2451515.501 


V 


.008 


38435 


0.00850 


0.00513 


BBSAG 121 


2451810.419 


V 


.005 


38855 


0.01103 


0.00762 


BBSAG 123 


2451841.3121 


cc 


.0058 


38899 


0.00823 


0.00481 


Zejda (2002) 


2451876.4201 


cc 


.0017 


38949 


0.00724 


0.00382 


Zejda (2002) 


2451878.533 


V 


.008 


38952 


0.01360 


0.01018 


BBSAG 124 


2451899.5900 


cc 


.0003 


38982 


0.00521 


0.00179 


Agerer & Hubscher (2002) 


2452190.295 


V 


.003 


39396 


0.00783 


0.00437 


BBSAG 126 


2452204.3365 


cc 


.0010 


39416 


0.00573 


0.00227 


BBSAG 126 


2452213.4637 


cc 


.0005 


39429 


0.00460 


0.00113 


BBSAG 127 


2452213.4662 


cc 


.0070 


39429 


0.00710 


0.00363 


Zejda (2004) 


2452260.516 


V 


.004 


39496 


0.01086 


0.00739 


BBSAG 127 


2452524.525 


V 


.007 


39872 


0.00030 


-0.00320 


Diethelm (2003) 


2452531.5507 


pe 


.0002 


39882 


0.00420 


0.00069 


Agerer & Hubscher (2003) 


2452885.446 


V 


.002 


40386 


0.00094 


-0.00261 


Diethelm (2004) 


2453302.5428 


ec 


.0010 


40980 


0.00301 


-0.00059 


Diethelm (2005) 


2453656.4422 


ec 


.0003 


41484 


0.00385 


0.00019 


Zejda et al. (2006) 


2453705.5937 


pe 


.0008 


41554 


0.00278 


-0.00088 


Hubscher et al. (2006) 


2453988.5713 


ec 


.0001 


41957 


0.00197 


-0.00172 


Brat et al. (2009) 


2454019.4675 


cc 


.0001 


42001 


0.00226 


-0.00144 


Brat et al. (2009) 


2455114.8667 


cc 


.0006 


43561 


0.00116 


-0.00268 


Diethelm (2010) 


2455149.2747 


cc 


.0003 


43610 


0.00236 


-0.00149 


this paper 


2455149.9768 


cc 


.0003 


43611 


0.00228 


-0.00157 


this paper 
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Table 4 Photometric Solutions for PS Per. 



^ ammeters 


r notometric eiements 


horrors 


9i = 92 


0.32 


assumed 


A! = A 2 


0.5 


assumed 


X\hol 


0.480 


assumed 


X2bol 


0.536 


assumed 


XXV 


0.506 


assumed 


XlR 


0.414 


assumed 


X2V 


0. IZO 


assumed 


X2R 


0.600 


assumed 


Tj. 


6750K 


assumed 


q (M 2 /Mi ) 


0.518 


0.003 


flin — 


2.8944 




Qout 


2.5907 




T 2 


4822K 


7K 


i 


89.°86 


0.°35 


n 2 

ri (pole) 


3.590 


0.008 


0.3229 


0.0008 


ri (side) 


0.3317 


0.0009 


r-i(back) 


0.3402 


0.0010 


r 2 (pole) 


0.3025 


0.0004 


r 2 (side) 


0.3159 


0.0005 


r 2 (back) 


0.3483 


0.0005 


E(O-C) 2 


0.0022 





